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A consideration of the lack of history dependence in the non-equilibrium steady state of a quantum 
system leads us to conjecture that in such a system, there is a set of quantum mechanical observables 
whose retarded response functions are insensitive to the arrow of time, and which consequently satisfy 
a quantum analog of the Onsager reciprocity relations. Systems which satisfy this conjecture can 
be described by an effective Free energy functional. We demonstrate that the conjecture holds in a 
resonant level model of a multi-lead quantum dot. 
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Although the fundamental principles of thermal equi- 
librium were established by Boltzmann more than a cen- 
tury ago, their generalization to the non-equilibrium 
steady state has proved elusive. The non-equilibrium 
steady state is thought to be denned by a set of character- 
istic variables such as the current, the thermal and chem- 
ical potential gradient and as such, it is expected to be 
independent of the history of how it was prepared. This 
has led to the notion that general principles should gov- 
ern the instantaneous properties of the steady state. One 
recurring idea is that a generalized free energy functional 
might apply to the non-equilibrium steady stateETQ. Thk 
was first-sneculated by Rayleigh in the late 19th century.EJ 
OnsagerB'0 later used his reciprocity relations to support 
this conjecture, but the idea has remained controversial 
to the present day. 

Non-equilibrium steady state behavior plays an impor- 
tant role in electronic transport theory, and becomes par- 
ticularly important in-driven nano devices, such as a d.c. 
biased quantum dot.Ll Variants on Rayleigh's approach 
would be invaluable in this new context, and might pro- 
vide an important first step along the road to ftaJtz- 
mann's approach the non-equilibrium steady state.ucl 

Recent work on non-equilibrium hydrodynamics has 
shown how Onsager's reciprocity relations .can be gener- 
alized to the non-equilibrium steady state.Hu This moti- 
vates us to re-examine Onsager's reciprocity relations in 
the context of non-equilibrium quantum physics. By con- 
sidering the history independence of the non-equilibrium 
steady state, we are led to conjecture that Onsager's 
reciprocity theorem continues within a limited class of 
quantum variables, in the non-equilibrium steady state. 
Within this restricted class of variables, the concept of a 
Free energy can be used to describe the steady state of 
non-equilibrium quantum systems. 

The lack of history dependence of the equilibrium 
steady state means that the work done on the system by 
coupling various internal degrees of freedom Ai (i = l,n) 
to corresponding external "forces" A, (t) , 



does not depend on the path P over which the (Ai) are 
adiabatically incremented to their final value. If we in- 
crement Xj(t) at two different times t 2 and t\ > t 2 , we 
may do it two ways, illustrated in Fig. (RT). 





FIG. 1. Two variations in the path P where the incre- 
ments in \j at times t% and t\ > ti are interchanged. 

In the first variation Aj(ti) — > \%{t\) + SXi and Xj(t 2 ) —* 
Xj(t 2 ) + SXj, whereas in the second the variations are 
reversed SXj «-> 5Xj. The second-order change in the 
work done along both paths must be equal, i.e 



S 2 W = 6X l SX i 



5X 3 (t 2 \ 

from which if follows that 



SXiSX. 



V SXi(t 2 ) 



5Xi(t 2 



8X 3 {t 2 ) 



(1) 



(2) 



We can relate these functional derivatives to the corre- 
sponding response functions, 



5{A 3 {t)) 



SXi(t') 

from which it follows that 



= -i{[A j (t),A i (t')])e(t-t') 



(3) 



W = J2 (Mt))dXi(t), 



- i ([A 3 (l),M2)])d(l 2) = Aj(2)])e(l - 2). (4) 

These are the quantum generalization of Onsager's reci- 
procity relations!™. The relations are understood to hold 
only in the long-time limit corresponding to a slow adi- 
abatic variation of the source terms. Onsager identified 
relations with the microscopic reversibility of the equa- 
tions of motion and the absence of any "arrow of time" 
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in thermal equilibrium. This derivation shows how reci- 
procity is directly related to a lack of history dependence. 
Since our proof makes no reference to thermal equilib- 
rium, it offers the intriguing prospect of an extension to 
the non-equilibrium steady state. 

To extend the discussion away from thermal equilib- 
rium, we considex-.a tiny system "S"' which may be a 
quantum dotu'EjO, a quantum wireEj, or other small 
system that is coupled to two very large baths of elec- 
trons ("leads") at different chemical potentials fiL and 
where Hl > /i_r. The entire coupled system is com- 
pletely isolated from the outside world. 
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FIG. 2. The non-equilibrium steady state is obtained 
by adiabatically connecting system S to two heat baths 
at chemical potentials ^l,r- 

If we connect S to the leads at time t = 0, then after an 
equilibration time n, the system will arrive at a steady 
state where a current flows from the left to the right- 
hand lead. This state persists for a long time T2(L) until 
a substantial fraction of the additional electrons on the 
left lead have flowed into the right lead. The time T2(L) 
will diverge rapidly as L — > oo, which permits us to define 
the steady state value of some variable A as 

(A) = lim (A(t)) 

L — >oo 

with the understanding that T2(L) >> t >> t\. 

Suppose the steady state is arrived by adiabatically 
turning on an interaction Hi = ghi between the leads, 
and by coupling source terms Xj to various quantities Aj 
which are localized within S. Since the combined system 
is closed, when we adiabatically change these variables 
the amount of work done in reaching the steady state is 
simply the change in the total energy of the system 

Wne = I \hi{t))dg{t) + (A^tydXi. 

If the work done Wne is independent of the path by 
which g and the Xj reach their final values, then we can 
use the previous proof to show that the corresponding 
variables satisfy a quantum reciprocity relation. The con- 
verse will also hold true. This motivates the "Quantum 
Reciprocity Conjecture" : 

In the non equilibrium steady state, the set 
of quantum mechanical observables contains 
a non-trivial subset V of "protected" quan- 
tum observables V — {ai, 02 ... , a„} whose 
correlation functions in the steady state are 
insensitive to the arrow of time, and which 
consequently satisfy a quantum mechanical 
analog of the Onsager reciprocity relations 

<[«(!), 6(2)]) = <[&(!), a(2)]>, (a, b e V). 



Of course we do not expect the reciprocity relation to 
extend to all variables, as it does in thermal equilibrium, 
because this would mean that the arrow of time is com- 
pletely invisible. 

Consider the retarded and advanced Green functions 
between protected variables, 



G 



(R,A) 



Ti([a(l),bp)])0±(t2-ti) 



(5) 



where 6±(t) = 8(±t). Since a and b are hermitian, these 
are real functions (G R ' A (t) = [G R > A (t)]*). The conjec- 
tured Onsager relations mean that in the steady state, 
they also satisfy 



G 



(R,A), 



tl) 



_ n (R,A) 



t 2 ), 



(t 2 -* 1 ) = G^(t 2 -t 1 ), 



(6) 



where the order of the subscripts and time variables is 
important. If we write G R (ti — £2) = [G R (ti — t2)] in 
the first relation, and then Fourier transform, we obtain 
the more familiar result 

G&M - G« 6 ( W )* 

which means that the retarded and advanced Green func- 
tions of protected variables share the same spectral de- 
composition 



G 



(R,A) 
ab 



dE 



1 



-E±i5 



Aab(E) 



where A ab (E) = ±lm[G^' R) (E)}. 

Provided that the set of protected quantum variables 
includes the interaction Hj = ghi, then we can define an 
effective Free energy from the virtual work done Wne 
in reaching the steady state. Suppose we evaluate Wne 
along the two paths shown in Fig. [| Since Wne is the 
same along both paths, for small A A we have 

f 92 da' 

A(g u X)AX+ / -^j-Hj(g', X + AX)dg 



f 92 da' 

A(g 2 ,X)AX+ / JL Hl (g',X)dg', 
Jgx 9 



so that 



where 



AA = A(g 2 ,X)-A(g 1 ,X) = — AF 



— f 92 % Hl {g',X) 

J 91 9' 



(7) 



(8) 



(9) 



Thus if reciprocity holds, the change in the variables {Aj} 
associated with a change in the coupling constant g can 
be evaluated as derivatives of a single Free energy variable 
AF. 

We now illustrate the correctness of this conjecture in 
a simple non-interacting model. We consider a single 
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X + Ak 







H(g , l + Al) 






'A(g,,X) 
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FIG. 3. Two paths for turning on the interaction and 
source terms. 

resonant level in a quantum dot carrying a D.C. current 
between two or more leads according, where the Hamil- 
tonian H = Hp, + Hj and 



E 



h = e ( k ) c Lk CT c ",k ff 

a,k<r tr 



Here a = l,iV labels the leads, each one characterized 
by a distinct chemical potential /x Q , tda — £d — &B is the 
energy of the localized state in the dot in a magnetic field 
B, J is the overall coupling constant and 7 Q is a parame- 
ter which sets the relative strength of hybridization with 
the a lead. This is_an exactly solvable problem, and has 
well known resultsEj found by the Keldysh method. 

As a first step, by comparing the retarded and ad- 
vanced correlation functions, we are able to explicity con- 
firm that the interaction, together with the dot magne- 
tization M and occupancy rid, form a set of protected 
variables {-£//, M, rid} which satisfy reciprocity and for 
which a Free energy functional can be defined. 

For example, to confirm the relation 

([H I (t 1 ),n(t 2 )]) = ([n(t 1 ),H J (t2)]), (10) 
we compare the retarded and advanced Green functions: 

de 



2tt 



[£ddt(e)(iri)£ Ccid t(e + w) 



and 

t-4 



S dr t (e)(iri)^ rfd t(e + w) 



Y~ [^c a dt(e + w)(in)^ dd t(e) 
+Sddt(e + ^)(""i)£ dc t (e) . 



The correlation function between M and rid identically 
vanishes, trivially satisfying reciprocity. 

We now confirm that an effective Free energy correctly 
determines the occupancies and magnetization. The ex- 
pectation value of the interaction energy determined by 
the equal time Keldysh Green functions between the con- 
duction and dot electron, given by 



where the G a b ■ refer to the Larkin-Ovchinikov matrix 
Greens functionliJIlj between electron fields a and b and 
the trace is over Keldysh indices. By writing these ex- 
pressions out explicity, we are able to explicitly confirm 
that they are related by complex conjugation, G^ in (uj) = 
[G'fj n (uj)]* , from which reciprocity between rid and Hj 
holds. A similar method enables us to check that 



[F J (ti),M(t 2 )]) = <[M(t 1 ),ff / (t 2 )]). 



(13) 



[G K 



After integrating over the coupling constant we obtain 

AF eff = 



J 

a, a 

+ 27rTlogr 



iA — ji a 



2niT 



(do 



2mT 



A In 



—) 

2ttTJ 



(14) 



where A = J2 a ^Pi^la) 2 - The expectation value of local 
state occupancy rid an d magnetization M are then 



i n d) 



dAF ( 



cir 



(M) = 



de d 

dAF cff 
dB 



Cl, 



C2, 



(15) 



where the constant terms gives the limiting value of the 
occupancy and magnetization when J — > 0. We can fix 
these constants by using the condition that (rid) — > 1 and 
(M) — > as A — > oo. which gives 



{n d ) = 1 + V ^Im 

7T 

a, a 

(M) = V ^alm 



Cdo - ^ a +iA 



1> 



2niT 
1 e d + oB - fl a 



iA 



2inT 



(16) 



(11) 



(12) 




FIG. 4. Distribution function of rid as a function of ej. 
/ii = 1, fj.2 = -1, Ai = 0.75, A 2 = 0.25, A = 0.01, and 
T = 0.001. 

Both results can be independently confirmed by direct 
calculation from the Keldysh Green functions. It is re- 
markable that the derivative of a single Free energy func- 
tional reproduces the results of two separate Keldysh cal- 
culations, even though a D.C. current is flowing through 



3 



the dot. It is interesting to see that even at the zero 
coupling limit, the occupancy and magnetization of the 
"dot" a non-thermalized form, and depends on the ra- 
tios between hybridization j a . The non-thermal function 
nd(ed) is reminiscent oil-thc occupancy observed in quan- 
tum wire experiments.E^I Here the parameters Xi play the 
similar role of distances between the measured point and 
leads in the experiment. 

It is instructive to examine the the magnetization in 
the two-lead case which for zero temperature is 



X(B,A) = 



2A{B 2 



V 2 



ir((B - V) 2 + A 2 )((B + V) 2 + A 2 ) 
whilst for A — > 0, 



X{B,T) 



1 

4T 



sech 



B + V 
2T 



sech 



B-V 
2T 



(17) 



(18) 



In both limits, the bias voltage dramatically reduces the 
susceptibility and at a finite voltage the T = mag- 
netic susceptibility in the limit of J — > is always zero. 
Non thermal magnetizations of this kind have recently 
obtained in the zero order -magnetic susceptibility calcu- 
lation for quantum dot£3ii3. Can we extend the set of 
"protected" variables to include other quantities of inter- 
est, such as the current or the spin current? The answer 
appears to be "no". When we directly compare the re- 
tarded and advanced correlators involving any operator 
that involves the lead electrons, other than Hi, we find 
that they are not complex conjugates. This means that 
we can not change the ratio of the couplings 7 a as we 
turn on the interaction, for to do this would be to in- 
troduce new variables which do not satisfy the Onsager 
reciprocity relation with hj. 

The validity of our conjecture in more complex systems 
is an open issue. We can not prove that reciprocity is sta- 
ble against the presence of interactions within the dot, 
but we have circumstantial support for this idea. The 
above methods can be used in the large N limit of the in- 
finite U Anderson model to examine how the mean-field 
equations evolve away from equilibrium. We have also 
compared the local susceptibility in the non-equilibrium 
Kondo problem obtained using the reciprocity conjecture 
with that obtained using Majorana techniqes.Ea An in- 
teresting recurring feature of these calculations, is the 
appearance of non-thermal distribution functions in the 
limit that the coupling with the leads is taken to zero. In 
interacting systems, these limiting distribution functions 
will need need to be computed self-consistently from the 
limiting form of the Dyson equation, before the change 
in Free energy can be computedE3 

In conclusion, we have examined the idea that the 
principle of virtual work can be extended to the non- 
equilibrium steady state of quantum systems. This has 
led us to conjecture the existence of a class of steady 
state variables which satisfy the quantum generalization 
of Onsager's reciprocity relation out of equilibrium. If 
this conjecture holds, then the notion of a free energy 



can be extended to the quantum non-equilibrium steady 
state, permitting the expectation values of steady state 
variables to be computed as derivatives of a free energy 
functional. This idea works for the simplest possible ex- 
ample, and leaves open the possibility that it will apply 
to more complex and interesting interacting situations. 
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